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We identify the independent dimension 6 twist 4 gluon operators and calculate their renormaliza¬ 
tion in the pure gauge theory. By constructing the renormalization group invariant combinations, 
we find the scale invariant condensates that can be estimated in nonperturbative calculations and 
used in QCD sum rules for heavy quark systems in medium. 
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I. INTRODUCTION 

Understanding the changes of the matrix elements of 
the gluon operators near the critical temperature in QCD 
offers a useful picture on the nature of the QCD phase 
transition[l|. These can also be used in QCD sum rule 
analysis to understand the changes and melting of heavy 
quark system at finite temperature [l-dj. 

In the pure gauge theory, the lowest dimensional op¬ 
erators are the scalar gluon condensate and the twist 
2 gluon operator. These dimension 4 operators can be 
reexpressed in terms of the electric condensate and the 
magnetic condensate. The temperature dependence of 
these operators can be calculated directly from lattice 
calculation of the space time and space space elementary 
plaquette[l], @ or from combining the calculation of the 
energy density and pressure The calculations show 
that while there is rapid change of the electric conden¬ 
sate across the phase transition temperature, the mag¬ 
netic condensate changes very little [Q, oj. 

Using the temperature dependence of the dimension 
4 condensates as the input in the QCD sum rule ap¬ 
proach for the heavy quark system, J/'ip and r] c was found 
to undergx) a rapid property change across the phase 
transition]^, Ej and to their dissociation |4j] slightly above 
the critical temperature. Moreover, it was also found 
that the free energy extracted from lattice calculation 
is the relevant potential to describe J/ip in a potential 
model[8j. The extension to finite density also has inter¬ 
esting application [9j. 

To further understand the phase transition in terms of 
local operators and to expand the findings for the char- 
monium system by using QCD sum rule to dimension 6 
level, we will identify the dimension 6 and twist 4 gluon 
operators and calculate their renormalization in the pure 
gauge theory. The renormalization of scalar dimension 
4 operators and scalar dimension 6 operators are well 
known (To], Hlj |. Our result completes the calculation of 
renormalization of all the dimension 6 gluon operator, 
hence will be a first step toward identifying their mixing 
and thus a systematic analysis in the operator product 
expansion (OPE) of heavy quark correlation functions up 
to dimension 6[12|. 

In section II, we will identify the independent operators 
at dimension 6. In section III, we will renormalize these 


independent operators up to one loop order. The scale 
invariant vacuum condensate will then be given in section 
IV. Section V is a summary. 

II. INDEPENDENT OPERATORS 

The gauge invariant dimension 6 operators are ob¬ 
tained by combining the covariant derivative D M and the 
field strength tensor G^. To find the independent even 
parity operators, we use the Bianchi identity and symme¬ 
try property of the indices. Here, we start from the oper¬ 
ators that are of the type ( D a Gbc){T>dG e f ); that is, mul¬ 
tiplication of two covariant component each composed of 
a covariant derivative acting on the field strength tensor. 

For the scalar operator, the indices 1 abodef' have to 
become l aabbcc ’ type. Considering the indices, the co¬ 
variant component ( DG ) can be a term in one of the 
two types of the Bianchi identities. In the first case, the 
three indices l abc ’ are independent while in the second 
case, two indices are identical and summed over l aab\ 

Type 1 : D a Gb c + DbG ca + D c G a b = 0 , 

Type 2 : D a G a b + DbG aa + D a Gb a = 0 . ( 1 ) 

The scalar dimension 6 operator can be obtained from 
one of the terms in the product of the same type in 
Eq. HI . Initially, four operators can be constructed from 
{Type 1) x {Type 1), and one operator from {Type 2) x 
{Type 2). However, among the four types of operators 
coming from {Type 1) x {Type 1), using the symmetry 
property of the indices, one can show that only one op¬ 
erator is independent, irrespective of the order of their 
indices. Therefore, there exists two independent scalar 
operators. That is, 

D a G%D a G a ^, D,G“ a D„G a va . (2) 

Using the equation of motion, the second operator can 
be written in terms of quark operator, which vanishes in 
the pure gauge theory. Using higher dimensional Bianchi 
identity of the form [D,[D,G]] = 0, one can show that 
the usually quoted scalar operator can be obtained by 
combining the two independent operators. 

rabcs~ia ^b ^c _ p. ^a j-. ^a _ f jn n C< a 

IJJ '~ r fiv'~ r /j.a Kjr va — '- t /ack '-'Va ^ . {O) 
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Similarly, for the spin 2 operators, the indices ‘ abode /’ 
become ‘ abccdd that is, ab indices remain free while cd 
indices are summed over. Then, there are four types of 
Bianchi identities that become relevant. 

Type 1 : D a Gbc + DbG ca + D c G a b = 0 
Type 2 : D a G c d + D c Gd a + DdG ac = 0 
Type 3 : D c Gdd + DdGdc + DdG c d = 0 
Type 4 : D a G cc + D c G ca + D c G ac = 0. (4) 

The full operator can be obtained from one of the 
terms in the multiplication of Type 1 to Type2> and 
Type2,Type4 to themselves. In this case, one operator 
is obtained from {Type 1) x {Type 3), two operators from 
{Type 2) x {Type 2), and one from {Type 4) x {Type 4). 
However, using the higher dimensional Bianchi identity, 
one can obtain a relation among the four operators. 
Hence, there are only three independent dimension 6 
spin2 gluon operators. In this work, we will use the fol¬ 
lowing set [13 : 

scalar : f abc G“„G b m Gl a , D^ a D,G a va 

spin2 : D a G'; tlJ D a G%; D V G^ DpG^D^. 

On the other hand, using the equation of motion, we 
find that only two gluon operators of dimension 6 remain 
in the pure gauge theory. These are f abc G“G l \ a G° a 
and D a G a ^D a G a ^. The latter operator is proportional 
to /“^G^G^G^ with two spin indices {a/3). Intro¬ 
ducing the color E and B fields, we find the off diagonal 
components are of the forms E^B b ± B^ or E^E^B^, the 
matrix elements of which vanish in the medium at rest 
due to rotational invariance. Therefore, the two inde¬ 
pendent dimension 6 operators in the pure gauge the¬ 
ory that remain and that constitute the diagonal compo¬ 
nents and the scalar operators are f abc B a ■ {B b x B c ) and 

jabc B a . ( E b x E cy 


To study the renormalization of the operator, we con¬ 
sider the following Green’s functions with external fields, 

<T“aUaOi> = Z 1A Z A {A;A b Al0 1B ) 

^ h ( 8 ) 
+ J2 z i,A^A b t/ A c x O jB ). 
j =2 

Here, is the background gluon held and Za the back- 
groundheld renormalization constant. Oib represents 
the bare operator with renormalized fields and coupling. 

The diagrams that contribute to the renormalization 
are shown in FIG.l with the Feynman rules given in 
FIG.2. The two gluon vertex comes from contraction 
with A c ^ l {p)A b l/ {q), the three gluon vertex comes from 
A/ l {p)A b {q)A , ^{r), and the four gluon vertex comes from 
A“{p)A b v {q)Al(r)Ai{k). 



FIG. 1: Diagrams contributing to the renormalization of Oi 
to one loop order in the pure gauge theory. 

The calculation is performed using dimensional regu¬ 
larization D = 4 — 2e in Feynman gauge for SU(N). The 
calculation for the other operators involves the same dia¬ 
grams and hence can be repeated similarly. The following 
is the collected result of the renormalization constants. 




III. RENORMALIZATION 

The renormalization of scalar operators are reported 
in Ref. 0. Here, we will focus on the spin2 trace¬ 
less (Twist4) part. We will use the three independent set 
as mentioned in the previous section after making the op¬ 
erators symmetric and traceless with respect to the two 
spin indices. We will therefore discuss the renormaliza¬ 
tion of the following three operators: 

Ol = DflGZrDcGZrlsT, (5) 

O 2 = D^GI^DMIAst, (6) 

0 3 = DpG^DM^sT, (7) 

where we have suppressed the external indices a, f3 in the 
left hand side and O a p\sT means that 1/2 {O a p + Op a ) — 
1 / 4g a pTr{0 a p). First, we will study the renormalization 
of 0\ up to one loop order using the background held 
method with zero momentum insertion fisj. 


Si.i = 1 + -- 1 


3 Nos 
4 Tie 


_ _N a s 
1,2 12 7T6 

2TV ct. 


■Zl,3 — 


■2^2,1 — 0 


Z2 ,2 — 1 


N_ Os 
3 ne 


%2,3 — 


N_Os_ 
24 7TC 


^3,1 — 0 

7 _ N a s 

^ 3,2 — - 

0 7T6 

Z -1 + 7iVas 
Z3 ’ 3 - 1+ Ti^- 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 

( 17 ) 
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4?“ pP S ab (p 2 g^ U - P^P U ) - g a ^ (4p 4 5 afe g /i,I/ — 4p 2 p^p U 5 ab ) 

4 



- ^af a bc( 2 P°‘P > '<l ,3 9 IJt ' ~ 2p a p l 'g^g X P + 2p^ q a (p X gP" - p" g X P) - p X g a 13 g P u p ■ g + p" g a P g X Pp ■ q + 2p a p X r? gP u + 4p a p X 
g^ U + 4 pP p X r^ g ai/ — 2p a p u g X ^ — 2 p X p u g a ^ + 2 p&r a (p X g tiV — p U g X ^) — p X g a ^ g^ p ■ r — 4p a g^ u g X ^ p ■ r 

— 4 p^ g aU g X ^ p ■ r + 3 p U g a ^ g X ^p • r) + (5 other terms from contraction order) 



- ~ 9 2 f a bxfcdx(. AkOL 9^^ 9 UU> + 4r a p@ g X ^ g UUJ +4 k a q^ g X ^ g UUJ +4r a q@ g X ^ g UUJ + 4k a g@^ q X g UUJ + 4r a g^^ q X g UUJ 
4 

+ 4g aX p^k+g^4g olX q^k^g I ' UJ - g*? p X g™ - 2 g a f* q X k» g™ - g a ^q X r^g^ - g^ g X » k ■ pg^ + 8g aX 

k ■ qg U “ - 3 g a Pg X f*k • qg™ - g° t ^g X ^p • rg vu - g a ^g X ^q ■ rg™ - 4 g aX p^g^k u - 4g aX q^ g^ k v + g a Pp X 
k u + g*? q X g»“ k v - 2g a ^p X g^p v + 2 g a P g X >*p" p“ + 4p a gP v {p X g» u - g X ^p UJ ) + 4g av p? (p A g»“ - g X »p“) 

- 4k a g^^ g Xu q^ — 4r a g&^ g Xu q u + g a & g Xu k^ q W + g a ^ g Xu q W — 8g aX k v q^ + 2g a @ g Xfl k u q^ ) + (23 other terms) 



. * 0 u X cxu u u X a/3 _ 2 A an 0u 0 2 A cxu 0u . 2 A a/3 ui/ , _ a u, 0u, A . A x . 0 Xu u 

— z 9f a bc( 2 P P ( l 9 —PP 9 9 — 2 p q g g’^ - 2p q g ^ g H + p q g H g^ +2 p p^ (g H ( q +r ) + r'g -r 

0X\ , A cku u, u\ cxB 0 8 u a A . A „ u, u <y.B <-,2 8 Xu <-,2 „ A olu Bu 0 2 A oiu Bu , 2 X 

g H ) + 2p H p^r g — p^ p r g H — 2p H p^r g + p p^ r g H — 2p r' g ^ g -2p r g g Hh ^ —2p r g ^ g H + p r 

cxB uu . 2 u cxB Xu . 0 2 u cxu 8X . 0 2 u cxX Bu 2 u cxB Xu . „ cx Xu, B u 2 Bu^ u cxB Xu 

9 9 +P ^ g 9 +2 p r g ^ g^ +2 p r g g ,h ^ - p rg H g^+2rg (p H p^ - p g H ^)~P^g H g p ■ r) 

+ (5 other terms) 



^ 2, , / a cx BX u uu , . cx X 8 u utu „ cxB X u uuj cxB, X u uu . cxu BX 2 uu . ct X Bu 

- 9 f abx f cdx ( 4 P 9 P 9 +4 g p H p^ g -2 g H p p^ g -g H k q^ g - 4g ^ g’ p g -4 g g H ^ 

4 

+ 2g°‘f 3 g X »p 2 g<'“ + g^g^q^k" + 4 g^k?g Xu p v - ig^g^k X p" + g a P k X g^p" - g a P g X “ p" + Ag a » k? 
g^g" - 4 g^g^^g" + g a f> k X g»“ - g a P g Xu k» q" - g^ q X g^ k“ + 9 “' 3 g X " q* fc" + 4 s «<‘/ i P I 'l.“ - g 

g X > X p‘ , k U + 4 S “'Wfc"’ - g al3 g XlL q V + g°‘f i k X g^q^ - g a P q X g^ r“ + g 01 ? g X “ q^ r“ + ig 1 *^ gP X p" r“ 

- g^g^p"^ + ig°‘»gl >X q ''r“ - g*? g X »q“ r“ + 4 q a gl" (k? g X “ - gP“ k X + g^ X (.h a + r U » - ig au q» (k? g X “ 

- gP^ k X + gP X (k^ + )) — g a ^ g XtAj g* LV k ■ q ) + (23 other terms) 


- s ab( 4 P a pP - p 2 a a ^)(p IJ 'p 1 ' - p 2 g^ u ) 



• c t a „<x Bu A u . n ocu /3 _ A u | . <xu B ~ A u . cx 8 Xu _u . cx0 „A u _u q0 A u_ u 0 cx8 _X u u 

— z 9 f abc ( 4r S' P P + 2 9 9 P + 4g r M p p^ + g H q p p^ + g H r p p^ - g 'p q p^ - 3g H p r p^ 

4 

„ otu 0X u „ ccA 0u u . 0 cx0 Xu u „ cxu 0 X u „ cxu 0 X u . 0 cxu 0 X u . cx 0u Xu 2 „ cxu 

— 2 g g^ p ■ rp^ — 2g g^ p ■ rp^ + 2g ^ g p ■ rp^ - 2g ^ p^ q p —2 g ^ p^ r p +2 g ^ p^ p r — 4r g^ g ^ p —2 g 

qPg X »p 2 - 4g otv g X ^ L p 2 - g^ q X g* v p 2 - g*? r X g^ p 2 + 2 g^g^r^p 2 + 2g aX g^ p 2 - 2g a ^g Xl 'r^p 2 + g*? g X » q" p 2 + 3 
g cx0Xp r u2 + 2q a g 0u (p X p p _ g Xp p 2 ) _ 2g cxp p 0 g Xu p . r + 2p “(2p^( g A ^^ + r X g^ 

+ g Xu - g X ^r u ) - g^ tJ, (q X p U + r X p u — p X r u + g Xl> p ■ r))) + (5 other terms) 



-9 2 fabxfcdx( 2 9 atJ 'P^ kX 9 l ' Ui + 4g atl q@ k X g uu + 2g a ^ p@ p X g"^ - g a>3 k X p^ g UUi - g a/3 p X p^ g UUJ - 2 g a 13 k X q^ g" u + g a/3 
4 

g X »p 2 g™ - 2 g^pPg^k" - 4 g a » g? g Xu k" + g a ^ g Xu k" + 2 g*? g X “ g* k» - 2g°‘» p? g X " k“ - 4 g^ q? g Xv k“ 

+ 2g°‘" gf 3 »q X k“ + 2g a »gP"q X k“ - g^ q X g^ fc“ + g“? g X “ fc" + 2®“? g Xv q» fc“ - 2g°‘ v g? X p»p“ - 2g aX g?" 

u u . cx 0 Xu u u „ cxu flu, A u „ an 0u, X u , a/3, A ui/ cj „ an 0 Xu u . cxu 0 Xu u . „ an flu 
V V +9 9 V P -2 g g^^k q - 2g ^ g H k q + g H k g^ q -2 g ^ p H g r -4 g ^ q H g r +2 g g H ^ 

X u , 0 cxu 0u X u a 0 X uu u . ext3 Xu u u , 0 ext3 Xu u u , . a flu,, A uu Xu , u Xu,,u , w XN 
9 r +2 g ^ g H q r - g H q g^ r +g H g p^r +2 g H g q^r + 4q g H ^(k g -g k -g (fc + r )) 

- 2p“(2p' 3 s A ' J 9 , '“ J + s' 3 ''(-fe^ g 1/ “ - p V“ + g X “k" + g Xv k“ + ")) + 2g au gP* g X “ k ■ q + 2g a » g?" g X “ k ■ q - g^ g X “ 

g^k.q + 2g“‘V X s '“*'p 2 + 2 s“ X s' 3 *'g^p 2 - s'** 3 g Xu g^ p 2 ) + (23 other terms) 


FIG. 2: Feynman rules in the background field method for the pure gauge theory. 
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IV. SCALE INVARIANT CONDENSATES 


The scale invariant condenstates can be obtained by 
diagonalization the following matrix Z. 


Z = 


i I 3Na s 
' 47T6 


0 

0 


Na s 


N, 

6ne 


3tt€ 


2 Na s 

i 


247T6 


(18) 


We then find the following new operator set, which cor¬ 
responds to the eigenvectors of Z. 


with corresponding factors of the coupling a s so that the 
renormalization of the coupling cancels that of the oper¬ 
ator [TH ]. 


(j)\ — (X s 11 (f^lnew) (25) 

4>2 = Ot s 1 2 3 4 5 6 ~^~ (0 2 new) (26) 

4>3 = &s 44 (03new) • (27) 


(Oinew) = (Oi> (19) 

{0 _ ) .(=551±|lA! 0l + i^ 0a + 0j ) (20) 
(o—> = ( - 653 + i± r g ° a +°») • 


These are renormalized multiplicatively without mixing. 
The renormalization constants correspond to the eigen¬ 
values of Z. 


(Oinew) — 
( 02 new) = 
(03 n ew) = 


i + ^)(o? newB ) 


47T£ 


1 + 


1 + 


(i5-Vrr)Na.\,„ 0 s 

48ne ) W 2newB / 

(15 + Vrf)Na s \ / ^ 0 \ 
487re I ■ 


( 22 ) 

(23) 

(24) 


V. SUMMARY 


We have identified and calculated the renormalization 
of the dimension 6 twist 4 gluon operators to one loop 
order in the pure gauge theory. Among the three inde¬ 
pendent operators, 0\ is related to the second moment 
of the usual dimension 4 gluon condensate and do not 
mix with other operators 0 2 and O 3 , which vanishes in 
the pure gauge theory. Hence, 0\ could be the first oper¬ 
ator that can be estimated in a non perturbative model 
or calculated on the lattice. With our calculation, the 
renormalization of all the dimension 6 operators are now 
known. The QCD sum rule methods for the heavy quark 
system in medium can now be systematically studied up 
to dimension 6 level. 
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